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In the momentum kick model, a near-side jet emerges near the surface, kicks medium partons, 
loses energy, and fragments into the trigger particle and fragmentation products. The kicked medium 
partons subsequently materialize as the observed ridge particles, which carry direct information on 
the magnitude of the momentum kick and the initial parton momentum distribution at the moment 
QQ , of jet- (medium parton) collisions. The initial parton momentum distribution extracted from the 

• STAR ridge data for central AuAu coUisions at y^ijvjv = 200 GeV has a thermal-like transverse 

' momentum distribution and a rapidity plateau structure with a relatively flat distribution at mid- 

, rapidity and sharp kinematic boundaries at large rapidities. Such a rapidity plateau structure may 

arise from particle production in flux tubes, as color charges and anti-color charges separate at 
high energies. The centrality dependence of the ridge yield and the degree of jet quenching can be 
consistently described by the momentum kick model. 
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I. INTRODUCTION 



In central high-energy heavy-ion collisions, jets are produced in nucleon-nucleon collisions and they interact with the 
dense medium produced in the interacting region. Depending on the relative azimuthal angle relative to the trigger 
pH ' particle, observed high-pt jets can be classified as near-side jets or away-side jets. An away-side jet is associated with 
a broad cone of particles pointing azimuthally opposite to the trigger particle direction. The strong attenuation of the 
away-side jet in its passage through the produced dense matter is one of the many notable experimental observations 
in relativistic heavy-ion collisions and is a signature for the production of the strongly-coupled quark-gluon plasma 

flSii. 

' On the other hand, a near-side jet is characterized by the presence of associated particles within a narrow azimuthal 
T-H ■ angle along the trigger particle direction. It retains many of the characteristics of the associated fragmentation 
' products as those of a jet in pp and peripheral heavy-ion collisions. The near-side jet occurs when the high-pt jet 
\^ , emerges near the surface of the producedparton medium. 

(3 ■ Recently, the STAR Collaboration d, |1, 0, H, H, [13, [HI, El, El El, [II observed a A0-Ar/ correlation of particles 
00 ! associated with a near-side, high-pt trigger particle in central AuAu collisions at y^SNN ~ 200 GeV at RHIC, where A(j) 
' and A77 are the azimuthal angle and pseudorapidity differences measured relative to the trigger particle, respectively. 
Particles associated with the near-side jet can be decomposed into a "jet component" , corresponding to fragmentation 
products of the near-side jet at {A(f>, Arj)^{0,Q), and a "ridge component" at Acfi^Q with a ridge structure in Arj. 
Similar A^-Ayy correlations associated with a near-side jet have also been observed by the PHENIX Collaboration 
^ ! and the PHOBOS Collaboration [l| . While many theoretical models [ll M, M, M, M, IH, M, M, 113, M, M, 

l30l. I31I. [3^. [33I. [33. [sH have been proposed to discuss the jet structure and related phenomena, the ridge phenomenon 
has not yet been fully understood. 

Previously, a momentum kick model was put forth to explain the ridge phenomenon [l3,|23,|2H. The model assumes 
that a near-side jet occurs near the surface, kicks medium partons, loses energy along its way, and fragments into 
the trigger and its associated fragmentation products (the "jet component") (Fig. 1). Those medium partons that 
are kicked by the jet acquire a momentum kick along the jet direction. They subsequently materialize by parton- 
hadron duality as ridge particles in the "ridge component" (Fig. 1). They carry direct information on the momentum 
distribution of the medium partons at the moment of jet-(medium parton) collisions, for which not much information 
has been obtained from direct experimental measurements. As the early state of the medium partons is an important 
physical quantity, it is therefore useful to examine the early parton momentum distribution using the momentum kick 
model. 

A previous momentum kick model analysis gave theoretical results in qualitative agreement with experimental data 
[T^ . We arrived at the interesting observation that at the moment of jet- (medium parton) collisions the parton 
transverse slope parameter T is slightly higher and the rapidity width substantially greater than corresponding 
quantities of their evolution products at the end-point of the nucleus-nucleus collision. We would like to refine the 
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model and give a quantitative comparison with experiment. We also wish to explore the early parton momentum 
distribution over a wider kinematic range of different transverse momenta and rapidities to search for interesting and 
novel features of the initial parton momentum distribution. 



Near-Side Jet Parton 




FIG. 1: (Color online) Schematic representation of the momentum kick model. A near-side jet parton (represented by adjoining 
thick arrows) occurs in a dense medium, whose partons are represented by solid circular points. The jet parton kicks many 
medium partons, loses energy, radiates, and fragments into the trigger particle and associated "jet component" particles. The 
medium partons that are kicked by the jet parton acquire a momentum kick along the jet direction and materialize as associated 
"ridge component" particles. Not shown in the figure is the away-side jet opposite to the near-side jet. 

We shall show that the extracted early parton distribution has a plateau rapidity structure. Rapidity distributions 
in the form of a plateau have been known in QCD particle production processes both experimentally and theoretically. 
Experimental evidence for a plateau rapidity distributions along the sphericity axis or the thrust axis has been found 
in TT^ production in e+-e~ annihilation [la, [13, [H, [H, [10, l4l| . The rapidity distributions for production also 
show a plateau structure with a depression in an extended region around y ^ 0. As a function of energy, the shape 
of the rapidity distribution for the sum of all particles produced in e+-e^ annihilation is a plateau with either a flat 
distribution or a small depression at y ~ [40| . The width of the rapidity plateau increases as a function of increasing 
center-of mass energy. 

Theoretically, the rapidity plateau structure has been known in many earlier investigations of particle production 
processes in QCD, when a quark pull s away from an anti-quark at high energies [i^, [Ij, [3, [H, [4^ . The theoretical 
basis in the work of Refs. [13, [isl. |44|. [ist comes from the approximate connection between QCD and QED2 [i^, [H, 
|49| . [soj . We would like to show more explicitly here how the transverse confinement in a flux tube allows one to 
establish an approximate connection between the field theories of QED2 and QCD in high energy processes. Using 
such a connection, we wish to review here how the rapidity plateau structure occurs when a color charge separates 
from an anti-color charge at high energies. 

The early stage of the nucleus-nucleus collision comprise of many simultaneous elementary particle production 
processes involving a quark pulling away from an anti-quark (or qq diquark) at high energies. As elementary processes 
lead to plateau rapidity distributions, the rapidity distribution of the medium partons at the early stage of the 
nucleus-nucleus collision should retain the rapidity plateau characteristics. 

The parton momentum distribution is only one aspect of the momentum kick model. The momentum loss of the 
jet in its passage through the medium is another important aspect. While many theoretical treatments of the jet 
quenching phenomenon have been presented previously [5l| . the investigation of the jet quenching phenomenon in 
connection with ridge particles associated with the jet will provide a different and complimentary perspective. The 
ridge yield and the quenching of the emerging jet will depend on the number of medium partons kicked by the jet 
along its way. A successful simultaneous description of the centrality dependence of the ridge yield and jet quenching 
will provide a consistent picture of the interaction between a jet and the medium. It will also pave the way for a 
future Monte Carlo implementation of the momentum kick model where many refinements and improvements can be 
included. 

In this paper, we shall limit our attention to particles associated with the near-side jet. In the context of the 
present investigation, the ridge particle momentum distributions in nucleus-nucleus collisions refer implicitly to those 
measured on a "per jet trigger" basis, unless indicated otherwise. For brevity of terminology, the term "jet" will be 
used both for the parent "jet parton" that passes through the medium and also for the daughter "jet component" of 
associated particles. The ambiguity of the meaning of the term can be easily resolved by context. 
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This paper is organized as follows. In Section II we summarize the momentum kick model and relate the ridge 
yield to the number of kicked medium partons. In Section III, wc give the relationship between the initial and final 
parton momentum distributions under the action of a momentum kick. In Section IV, we specify how the initial 
parton momentum distribution is parametrized. In order to determine the initial parton momentum distribution from 
the observed total particle distribution in central AuAu collisions, we parametrize in Section V the jet momentum 
distribution associated with the near-side jet in pp collisions. The momentum distribution of mid-rapidity associated 
particles in central AuAu collisions at ^/snn = 200 GeV are then analyzed in Section VI. In Section VII, we display 
explicitly the initial parton momentum distribution extracted from the ridge data. In Section VIII, we examine the 
connection between QCD and QED2 in the presence of transverse confinement and study the origin of the rapidity 
plateau when a color charge separates from an anti-color charge at high energies. In Section IX, the field theory of 
bosonized QED2 is then used to study particle production as an initial- value problem. The evolution of the medium 
parton momentum distribution is discussed in Section X. In Section XI, we turn our attention to the propagation of 
the jet and the dependence of the jet fragmentation function on the number of jet-(medium parton) collisions. In 
Section XII, the centrality dependence the ridge yield and jet quenching is examined. In Section XIII, we examine 
the dependence of the ridge yield on the energy and mass number of the colliding nuclei. In Section XIV, we present 
our conclusions and discussions. 



II. THE MOMENTUM KICK MODEL 



It should be pointed out on the outset that the interaction between a jet and the medium can be described by 
representing the medium either as fields or as particles. In our momentum kick model, we choose to represent the 
medium as particles. We describe the interaction between the jet and the medium in terms of jet-(medium parton) 
collisions, from which each collided medium parton receives a momentum kick and subsequently materializes as a 
ridge particle. We have been guided to such a particle description because of the strong color screening in a dense 
color medium [s^. [ssj. The presence of the azimuthal kinematic correlation between the jet and the ridge particles 
lends additional support to the concept of jet-(medium parton) collisions as a central element of the phenomenon. 

As depicted in Fig. 1, the main contents of the momentum kick model can be briefly summarized as follows: 

1. A near-side jet parton emerges near the medium surface and the jet parton collides with medium partons on its 
way to the detector. It loses energy by collisions and gluon radiation. It subsequently fragments into the trigger 
particle and other associated fragmentation products. 

2. Each jet-(medium parton) collision imparts a momentum kick q to the initial medium parton of momentum pi 
in the general direction of the trigger particle to change it to the final parton momentum p = p/ = Pi + q, and 
it modifies the normalized initial parton momentum distribution dF/dpi to become the final parton momentum 
distribution dF/dpf. The kicked partons subsequently materialize by parton-hadron duality as ridge particles. 

Based on the above picture, we can describe the jet and the kicked particles in quantitative terms. We consider 
a nucleus-nucleus collision at a given impact parameter b with iVbin binary nucleon-nucleon collisions, and we label 
a binary collision by the index i. For the ith binary collision, there is a jet parton distribution dN^/dpj where the 
subscript j stands for the "jet parton" . The sum over all binary collisions for a given impact parameter leads to the 
total jet parton distribution dNj/dpj defined by 

dN^^^dNl 

dPj jr{ dp J ■ 

In a single pp collision, the yield of a trigger particle with momentum Ptrig is 

dNPV f dNPP- 

-T- (Ptrig)= / dpj——D{ptrig,Pj), (2) 

"Ptrig J aPj 

where D(ptrig, Pj) is the fragmentation function for fragmenting a trigger hadron of momentum ptrig out of a parent 
jet parton of momentum pj. For convenience of accounting in nucleus- nucleus collisions for a fixed ptrig, we rescale 
the fragmentation function by dividing the above equation by the quantity on the lefthand side, [c?A^^^/dptrig(ptrig)], 
to change the above equation to 



dNPP 

dpj-r—D{ptrig,Pj), (3) 
apj 
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where the re-normahzed fragmentation function -D(ptrig, Pj) is 

^(Ptdg, P,) = ^(Ptrig, Pj)/M^''7dPtrig(Ptrig)]. (4) 

Using the fragmentation function normahzed in this manner, a binary nucleon-nucleon colhsion (the z-th binary 
colhsion, say) produces a single trigger particle at the momentum ptrig, in the absence of jet-medium interactions, 

f dm 

/ dpj—-^D{ptiig;Pj) = (unquenched A'^trig arising from the ith binary colhsion) = 1. (5) 
J "Pj 

With the total jet source distribution dNj/dpj coming from all binary collisions in a nucleus- nucleus collision, we 
have 

Nbin= [ dpj^^D{puigPj)- (6) 



In the presence of jet-medium interactions, the total number of trigger particles iVtiig with momentum ptrig is 

^trig= / dp,^ J2 ^(^)e"'^"^i?(Ptrig; Pj^J2 - ^-)' (7) 

N=a n=l 

where N is the number of medium partons kicked by a jet of momentum pj along its way, A'max is the maximum 
N considered, and P{N) is a probability distribution of N, normalized by X)jv='o ^(-^) = 1- '^^'^ factor e^''"^ 
describes the absorption of the jet due to the inelastic fraction of jet-(medium parton) collisions. The quantity q„ is 
the momentum kick on the medium parton due to the nth jet- (medium parton) collision, and is the jet momentum 
loss due to the gluon radiation of the jet. We shall postpone our discussion of P{N) to Section XII. It suffices to 
indicate here that P{N) depends on the medium parton density along the trajectory and the jet-(medium parton) 
scattering cross section. 

As is obvious from Eq. ([7]), the number of trigger particles A^trig in a nucleus- nucleus collision (with the momentum 
Ptrig) will be equal to the number of binary collisions A^bin in the absence of any jet-medium interaction, 

A^trig ({Ca, q«, = 0) = TVbin. (8) 

The ratio of 7Vtrig({Ca, Qn, Ar} 7^ 0) in a nucleus- nucleus colhsion with respect to A'^trig({Ca, Qn, A^} = 0) in the 
absence of any jet-medium interaction is the Raa measure of jet quenching, 

. . r7/V ^ 
Raa = J^ dp,-^ J2 i^We-'^"^I?(ptHg;p,-Eq" ^r). (9) 

bmj Pj 

Because the kicked partons are identified as ridge particles by parton-hadron duality and two-third of the produced 
hadrons are charged, the distribution of associated ridge particles is therefore 



"-'^ ridge 

dp 



- I rfp,^ E''^^We-^"^^(Ptrig;P,-f:q„-A.)|^f]/«„^ 



(10) 



where < fun < 1 is the ridge attenuation factor for the n-th kicked parton to reach the detector and dFn/dp is 
the normalized momentum distribution of the n-th kicked medium parton, normalized to J dpdFn/dp = 1. We note 
that the righthand sides of Eq. ([7]) and (|10p differ only by the quantity in the curly bracket. It is convenient to 
define the expectation value (O) of a quantity O in the presence of the jet distribution, jet momentum loss, and jet 
fragmentation by 



{O) = I dp,^ ^ F(iV)e-';"^I?(ptrig; P, - q" - ^r)0 

N=l n=l 
- / E ^(A^)e-«"^i^(Ptrig; P, - E q" - ^r). (11) 

AT 1 ^ 1 
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Using this definition, the momentum distribution of the ridge particle momentum distribution per trigger particle is 
then the expectation value of the sum of the final momentum distribution of the kicked medium partons: 

TVtng dp \3^^ dp / 
The above equation can be re-written as 

{.fR){N)& + (j^ifRn - (/fi))[^(q„) - (^)]\| , (13) 



iVtrig dp 3 1' dp \ ^ ' dp 

\n— 1 , y 

where fj^ is the average attenuation factor, fj^ — fnRjN , and {N) is the expectation value of the total number 

of kicked medium partons per trigger particle. The quantity {N) is also the expectation value of the number of 
jet-(medium parton) collisions per trigger particle. We shall often label {N) alternatively as {Nk) with the subscript 
k to emphasize that this is the averaged number of fcicked medium partons per trigger. 
As defined by Eq. (fTTj) . {N) and {Nk) are given by 

1 r ^ 

{N) EE (Nk) dp.-^Y. NP{N)e-<^^D{pu,,; p, - q„ - A.). (14) 

trig J Pj ^j^J^ 

To understand the gross features of the phenomenon, we neglect the second term in the curly bracket of Eq. (|13p 
which arises from the fiuctuation of the quantities from their mean values. The formulation can also be simplified by 
taking the different momentum kicks q„ to be the average q. Using these simplifying assumptions, we then obtain 

' . (/«)|(A.0(^?)^ (15) 



A^trig dp 3 dp' 

Thus, the ridge particle distribution is separated into a geometry-dependent part (/j^)(2/3)(Afc) and the average 
normalized momentum distribution of ridge particles, (dF/dp). For brevity of notation, the bracket symbol, (), 
for (dF/dp) will be made implicit, and the normalized ridge momentum distribution dF/dp will be understood to 
represent the average over the jet source distribution, jet collision locations, and jet energies. 

If one is interested in the total ridge yield by integrating over the ridge particle momentum, we then get 

^-{fH)UNk). (16) 

iVtrig O 

Our strategy is to study first the case of the most-central AuAu collisions at y^SNN = 200 GeV where the momentum 
distribution of the ridge particles and the average number of kicked medium partons can be inferred from experimental 
data [1, 0, 0- In Sections XI, XII, and XIII, we shall then examine the average number of kicked medium partons 
and the experimental ridge yield as a function of centrality, collision energies, and nuclear mass numbers (tI. [sI. [l^ . 
using the number of kicked medium partons for the most-central AuAu collision at ^^snn ~ 200 GeV as a reference. 



III. RELATION BETWEEN THE INITIAL AND FINAL MOMENTUM DISTRIBUTIONS 

In the momentum kick model, the normalized final parton momentum distribution EdF/dp at p is related to the 
normalized initial parton momentum distribution EidF/dpi at p^ at a shifted momentum, Pi = p — q, and we have 

m 



dF 



Ptdptdrjdcj) 



dF 



E 



Ptidpudyid(t)i Ei 



'1- 



{m? + p1) cosh y 



(17) 



where the factor E / Ei insures conservation of particle numbers and the last factor changes the rapidity distribution 
to the pseudorapidity distribution [4^ . The momentum kick q is expected to lie within a narrow cone in the trigger 
particle direction for a high-energy jet. To minimize the number of parameters, we approximate q to lie along the 
trigger particle direction. 
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To relate the final parton momentum distribution to the observed hadron momentum distribution, we assume 
hadron-parton duality which is a reasonable description for the hadronization of energetic partons. The final parton 
momentum distribution Eq. pT|) . multiplied by {fji){2/3){Nk), can then be identified with the observed (charged) 
hadron associated particle momentum distribution per trigger, dNch/Ntngdr]d(f>ptdpt, as given by Eq. (|15p . By a simple 
change of variables, we can further obtain dNch/ NtrigdArjdAcjjptdpt in terms of Arj = — ?/jct and A(j) = (j)—(j)jet-, relative 
to the trigger particle. The basic ingredients of the momentum kick model are then the magnitude of the momentum 
kick q, the normalized initial parton momentum distribution dF/dpi, and the average number of jet- (medium parton) 
collisions per jet. For numerical calculations, we set m = TO^r- 

The initial and final parton momenta can be represented in terms of Cartesian components in the collider frame, 
p ~ (pi,P2,P3), with a longitudinal p^ component, a transverse pi component, and another transverse p2 component 
perpendicular to both pi and p^. The coordinate axes can be so chosen that the trigger jet lies in the pi-ps plane. 
The initial parton momentum p,j = (j>ii,Pi2TPa) is related to the final momentum py = (p/i,P/2,P/3) and the trigger 
jet rapidity Tjjet by 

Pii = Pfi 

Pi2 = Pf2, 

Pa = Pf3 

For a given trigger particle pseudorapidity, these relations allow one to obtain p^ from p/ = p for the evaluation of 
the ridge yield per trigger particle. 



cosh ?7jet 

qsinh77jot 



cosh 



(18a) 
(18b) 
(18c) 



IV. PARAMETRIZATION OF THE INITIAL PARTON MOMENTUM DISTRIBUTION 

As the jet-(medium parton) collisions take place at different spatial and temporal locations during the passage of 
the near-side jet through the medium, the initial momentum distribution EidF/dpi in Eq. p7|) refers actually to 
an average over spatial and temporal regions during the early stage of the nucleus-nucleus collision. The 'initial' 
momentum distribution can also be called the 'early' parton momentum distribution. This initial parton momentum 
distribution EidF/dpi of the medium partons at the time of jet-(medium parton) collisions is not yet a quantity that 
can be obtained from first principles of QCD, although some of its qualitative features can be inferred from basic 
principles as will be discussed in Sections VIII and IX. Furthermore, jets occur at an early stage of the nucleus-nucleus 
collisions, whereas the momentum distribution of the bulk medium pertains to the state of the medium at the end- 
point of the nucleus- nucleus collision. Therefore the early parton momentum distribution near the beginning stage of 
the nucleus- nucleus collision needs not be the same as that of the bulk matter at the end-point of the nucleus- nucleus 
collision. 

Under the circumstances, the parton momentum distribution at the early stage of the nucleus-nucleus collision 
can only be obtained phenomenologically from the ridge particle data by representing it as a simple function whose 
distinct characteristics can be determined by comparison with experimental ridge data. 

The initial momentum distribution was parametrized previously as e~^' ^^'^h expj— -^/m^ -I- p^jT}/ y/w? + with 
m taken to be ttItt [19| . Although this is adequate for mid-rapidity and high-pf ridge particles [l3 |. it leads to too 
large a ridge yield both at pt ~ 1 GeV (dotted curve in Fig. 2) and at forward rapidities. Our understanding of 
the behavior of the early parton transverse momentum distribution has not reached such a stage that we can predict 
its low pt behavior definitively. If the partons arise from a deconfined medium with a finite transverse boundary, 
then the transverse parton momentum distribution at small pt will be fiattened from an exponential distribution, as 
shown in Figs. 1 and 2 of [53 |. Transverse distributions of this type can be described by replacing the denominator 
^/Tr^?+]^^ with where m^j can be adjusted to lead to the correct ridge yield aX pt ^ 1 GeV. The extracted 

transverse momentum distribution may provide useful information to study the transverse radius of the deconfined 
parton medium [H, [IHl ■ 

The difficulty with the forward rapidity region can be resolved by noting that the Gaussian rapidity distribution of 
[l^ does not take into account the kinematic boundary restrictions on phase space. The large values of ay extracted 
from the mid-rapidity data in [l^ imply that the rapidity distribution is quite flat in the mid-rapidity region. We 
can use a rapidity distribution that retains the flatness at mid-rapidity but also respects the kinematic boundaries at 
large rapidities and large pt ■ Accordingly, we parametrize the normalized initial parton momentum distribution as 

dF e-V"'+P?./T 
, , , . = ^■idgc(l~x)° . , . (19) 

pudptidy^d(j)^ y/mj + pj^ 
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FIG. 2: (Color online) The symbols represent STAR experimental data 0] and the curves theoretical results of 
dNch/NtiigPtdpt, for pp and central AuAu collisions. 



where A^dge is a normalization constant defined (and determined numerically) by 

dyid4>iPudpuAMgc{^ - xf 
In Eq. (fT9|) . X is the light-cone variable (isj l 



^,exp{-^^^^/T} 



(20) 



(21) 



rrib 

a is the fall-off parameter that specifies the rate of decrease of the distribution as x approaches unity, j/f, is the beam 
parton rapidity, and rrib is the mass of the bearn parton whose collision and separation lead to the inside-outside 
cascade picture of particle production [H, H, 0, 13 . As x < 1, there is a kinematic boundary that is a function 

of j/j and pti, 



^Vb-\yi\ 



(22) 



We expect to have a distribution centered around the nucleon rapidity, y^v — cosh^^(y^¥^/2mAr). For lack 
of a definitive determination, we shall set j/f, equal to and equal to ttIt^, pending their future experimental 
determination by examining the ridge boundaries. This form of the initial parton distribution leads to a restricted 
phase space that is smaller than that for a Gaussian rapidity distribution. As a consequence, it can lead to a smaller 
associated particle yield that agrees with experimental forward rapidity data as shown in Section VI. 



PARTICLE MOMENTUM DISTRIBUTION OF THE JET COMPONENT 



As a jet passes through the parton medium, the medium partons kicked by the jet will materialize to become 
particles in the associated "ridge component" , while the jet will fragment and radiate into the trigger and associated 
"jet component" particles. If the contribution from the jet component is known, we can separate out the ridge 
component using experimental data of total associated particles. The distribution of the "jet component" of (charged) 
associated fragmentation products is given by 



dp 



dN, 



P(A^)e-^"^i?2(Ptrig,P;P, 



N 

n=l 



A.), 



(23) 



where D2{pa, Pb! Pc) is the double fragmentation function for fragmenting into hadrons of momentum and p5 from 
a jet parton of momentum Pc. Fragmentation measurements Q suggest an approximate scaling relation 



-D2(Ptrig,P;Pc) ~ £'(ptrig;Pc)£)^(P;Ptrig), 



(24) 
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where i'zlPSPtrig) is approximately the same (within a factor of about 0.6 to 1.2) for dAu and AuAu coUisions in 
2.5 < pt.trig < 6 GeV (Fig. 5b and 5c of By applying this approximate scaling relation to Eq. (|^ and using 
Eq. ([7]), the (charged) jet component in an AA central collision per trigger is 



1 



Nti-ig dp 



^z(P;Ptrig) 



dp 



(25) 



Because of the approximate nature of the above relation (up to a factor of about 0.6 to 1.2), we need to make a 
quantitative check. In the region where the jet component has a prominent appearance, as in Fig. 3(d) for 2 < < 4 
GeV, the total dNch/NtrigdArj distribution at A77 ~ has indeed a shape similar to, but a peak magnitude about 
equal to, the pp near-side jet distribution. The total yield is the sum of the jet component and the ridge yield, and the 
ridge yield at A77 ~ is non-zero and nearly flat (Fig. 3(d)). The near-side jet component in AuAu central collisions 
per trigger is thus an attenuated pp near-side jet distribution, as expected for production in an interacting medium. 
If one assumes that fragmentation products lying deeper than an absorption length from the surface are all absorbed, 
then the average jet fragment attenuation factor is fj = e^^^^dx/ X = 0.632, which leads semi-empirically to a 
reasonable description of the experimental data as indicated below in Figs. 2 and 3. 
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FIG. 3: (Color online) The symbols represent experimental data and the curves theoretical results, for pp and central AuAu 
collisions, (a) and (b) give the dNch/Ntvi^dAtl) distributions, (c) and (d) give the dNch/Ntii^dArj distributions. 



The sum of the distributions (jl5p and (|25p . relative to the trigger particle rjjct and (pjct, is therefore given more 
precisely as 



TVtrig ptdptdAr]dA(j} 



AA 



total 



2 dF 



AA 



ridge 



fj 



ptdpt dArj dA(j) 



AA 



JCt 



(26) 



The experimental momentum distribution of (charged) near-side particles associated with the trigger in a pp colli- 
sion, measured relative to the trigger jet, can be parametrized as 



Ptdpt dArj dA(j) 



cxp{(m~ v/m2+p^)/rjct} 1 

Tjet(m + Tjet) 2TTaj 



-[(A0)2-H(AT,)^]/2a2 



(27) 



where A'jot is the number of (charged) near-side jet particles in a pp collision, and Tjct is the jet inverse slope 
("temperature") parameter. The above functional form of the jet fragmentation product cone in terms of Ai^ and A77 
was chosen because p2 = Pt sin A(j) and = pt sinh A77, and the square of the momentum perpendicular to the jet 
direction has a magnitude 



pI+pI= pI sin^ Acj) + pi sinh^ Aij ~ pl[{A(j)Y + {Arjf] for smaU Acj) and Aij. 



(28) 
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The above equation indicates the symmetry between A(/) and A?/ for a narrow jet cone. In this functional form of Eq. 
(P7)) for the jet cone, the width in A77 is equal to the width in A<j>. 
In our search for parameter values we find that the width parameter cr^ depends slightly on pt which we parametrize 

as 

(T4,=a^o I ° ■ (29) 

Experimental data for near-side jet particles inpp and central AuAu collisions obtained by the STAR Collaboration, 
within the detector acceptance of Ij^associatcdl < 1 and |77jet| < 0.7, are given in Figs. 2 and 3 0]. Figure 2 gives the 
dNch_/Nt^igPtdpt data, obtained by integrating dNch/NtrigPtdptdA(l)dAr] over the domain of \Ari\ < 1.4 and |A0| < 1.0. 
Figures 3(a) and 3(b) give dNch/ Nt^igdAcj) data, and Figures 3(c) and 3(d) give dNch/NnigdAr] data. They are obtained 
by integrating dNch/N^i-igPtdptdAc/jdAr] over the domains indicated in the figures. Specifically, Fig. 3(a) covers the 
domain of |?/| < 1 and 0.15 < < 4 GeV, 3(b) the domain of |?/| < 1 and 2 < pt < i, 3(c) the domain of \A(j>\ < 0.5 
and 0.15 < Pf < 4 GeV, and finally 3(d) the domain of \ A(j)\ < 0.5 and 2 < pt < 4 GeV. The domains of integration in 
a pp collision and a nucleus-nucleus collision are the same, and the distribution in A77 has been corrected for detector 
acceptance. 

The set of experimental pp near-side jet data of dN^^/ptdpt, dN^^ /dAcj), and dN^^/dArj, represented by open circle 
points in Figs. 2 and 3, can be described by Eqs. ([27|) and ([29|) . with the following parameters 

Tjet = 0.55 GeV, a^o = 0.50, = 1.1 GeV, and iVjet = 0.75. (30) 

Theoretical pp jet results obtained with this set of parameters within the specified experimental domain are shown as 
the dash-dot curves in Figs. 2 and 3. They yield a reasonable description of the experimental momentum distributions 
of jet particles associated with the near-side jet in a pp collision. 



VI. COMPARISON OF THEORETICAL NEAR-SIDE ASSOCIATED PARTICLE YIELDS WITH 

EXPERIMENT 

Theoretical evaluation of both the jet component and the ridge component for central AuAu collisions allows one 
to determine the total yield of associated particles as determined by Eq. ([26]) . A self-consistent search for the initial 
parton momentum distribution in Eqs. (|19p and (jl7p can be made by comparing the momentum kick model results 
with experimental data for dNch/NtrigPtdpt, dNch/NtrigdA(p, and dNch/ Ntj-igdArj for mid-rapidities in Figs. 2 and 3, 
and dNch/NtrigdA(f> for forward rapidities in Fig. 4. We find that the totality of the STAR associated particle data, 
from Pt = 0.15 GeV to 4 GeV and ?/ from zero up to 3.9 in central AuAu collisions at ^snn = 200 GeV [1,0,0, can 
be described by Eqs. (|26p and ([TT]) with parameters 

<Z = 1.0GeV, and (/i?) (iVfe) = 3.8, (31) 

in conjunction with the initial parton momentum distribution Eq. p9p with parameters 

T = 0.50 GeV, ma = 1 GeV, and a = 0.5. (32) 

We now discuss the comparison of the experimental data with theoretical results involved in the above analysis. 
In Fig. 2 the STAR experimental dNch/NuigPtdpt data are represented by solid circular points for central 
AuAu collisions and by open circle points for pp collisions. The theoretical results for pp collisions obtained with 
the parametrization of Eq. p7|) with parameters in Eq. ((30l) are shown as the dash-dot curves, which agree with 
pp near-side data in all experimental kinematic regions. Experimental ridge dNch / NtrigPtdpt data in central AuAu 
collisions 0] are also shown as solid squares and they are calibrated by using the data of Fig. 2 of 01- The solid curve 
is the theoretical result for dNch/NtrigPtdpt for central AuAu collisions, as the sum of the jet part and the ridge part, 
with the ridge part of the contribution shown as the dashed curve. They have been calculated with = 1 GeV. If 
we set md equal to to^, then we will get the ridge yield represented by the dotted curve, which over- predicts the ridge 
yield at ~ 1 GeV. 

Fig. 2 shows good agreement between theoretical dNch/NtngPtdpt with experimental data for central AuAu col- 
lisions. The theoretical transverse momentum distribution of the jet and the ridge components have very different 
shapes in the low pt region. The jet component dNch/Nn-igPtdpt decreases exponentially as a function of increasing pt- 
On the other hand, the magnitude of the final transverse parton momentum ptj is greater than the initial transverse 
momentum pu approximately by the amount q. The initial momentum distribution dF/pudpu has a peak at pu = 0. 
As a consequence, the theoretical ridge yield of final partons, dNf /NtngPtfdptf, has a peak around ptf q ^ 1 GeV 
and it decreases significantly for small values of pt, in contrast to the exponential behavior of the jet component. 
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It is interesting to note that the theoretical ratio of the jet yield to the ridge yield is greater than 1 for pt < 0.6 GeV, 
but is less than 1 in the interval 0.6 ^ Pt ^ 3.7 GeV. The ratio reverts to become greater than 1 at 3.7 GeV < pt- 
The change of the dominance of the ridge component as pt changes may lead to experimentally observable variations 
of the shape of the dNch/NtiigdArj as a function of pt- 
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FIG. 4: (Color online) Azimuthal angular distribution data at forward pseudorapidities for central AuAu collisions from the 
STAR Collaboration [0], compared with theoretical results shown as solid curves from the momentum kick model, (a) is for 
0.20 < Pt < 2 GeV, and (b) is for 1 < Pi < 3 GeV. 

In Fig. 3, the experimental total associated particle yields [1,0 are represented by sohd circular points for central 
AuAu collisions and by open circles for pp collisions. The theoretical results for pp collisions obtained with the 
parametrization of Eq. (j27p are shown as the dash-dot curves which agree with experimental pp near-side data. In 
these figures, the theoretical total yield and the ridge yield for central AuAu collisions are represented by solid and 
dashed curves, respectively. Comparison of the theoretical total yield and the experimental total associated particle 
yield for central AuAu collisions indicates general agreement over all azimuthal angles [Figs. 3(a) and 3(b)] and over 
all pseudorapidities [Figs. 3(c) and 3(d)], for both 0.15 < pt < 4 GeV [Figs. 3(a) and 3(c)] and 2 < pt < 4 GeV [Figs. 
3(b) and 3(d)]. 

One notes from Fig. 3(a) that for the region of 0.15 < pt < 4 GeV, which receives the dominant contributions 
from the low pt region, the widths of the azimuthal angular distributions of the ridge and the jet components are 
nearly the same, with the magnitude of the ridge yield slightly higher than the pp yield. On the other hand, in the 
region of 2 < pt < 4 GeV in Fig. 3(b), which receives the dominant contributions from the region near p^ '--^ 2 GeV, 
the azimuthal angular distributions of the jet component is narrower than the ridge component azimuthal angular 
distribution. 

We observe from Fig. 3(c) that the theoretical AuAu jet and ridge components have different shapes in 
dNch/NtrigdArj. The jet component maintains a sharp peak in dNch/NtiigdAr]. In the low pt region, the pseu- 
dorapidity distribution of the theoretical ridge component is significantly broader than the jet component and its 
magnitude remains to have a non-zero value at large \ Ar]\. In the high pt region in Fig. 3(d), the ridge pseudorapidity 
distribution is essentially flat and non-zero. The broad peak structure for the low pt region comes from the factor 
E/Ei in Eq. (fT7|) . arising from the difference of the momenta of the parton before and after the collision. This factor 
is close to 1 for the high pt region, and the flatness of the distribution is a reflection of the initial rapidity distribution. 

We turn now to forward rapidities where preliminary experimental data for central AuAu collisions have been 
obtained for 2.7 < [t^I < 3.9 (9|. We note that dV^^ch/A^trigdA^dAr] at A(/> ~ for \r]\ < 1 in Fig. 3(a) is an order of 
magnitude greater than the corresponding dNch/ NtrigdA(j)dAri for 2.7<|77|<3.9 in Fig. 4(a). This implies a substantial 
fall-off of the ridge yield dNch/NtYigA<f>dAri at A0~O in going from mid-rapidities to large rapidities. Even though the 
mid-rapidity data place a constraint on the flatness of the mid-rapidity distribution, they do not otherwise constrain 
the rate of fall-off of the distribution in the forward rapidity region. Measurements at forward rapidities in Fig. 4 
contain events with large 77 and p* that are either already outside the kinematic limits or close to the kinematic limits. 
Therefore, even with large errors, the forward rapidity data in Fig. 4 are sensitive to the constraint of the kinematic 
limits and the rate of fall-off of the initial parton momentum distribution as specified by the fall-off parameter a. The 
data of Fig. 4 lead to a = 0.5. 
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FIG. 5: (Color online) The momentum distribution of the associated particles as a function of the pseudorapidity relative to 
the jet pseudorapidity A77. The solid circular points are the data from the PHOBOS Collaboration [l^ and the curves give 
theoretical predictions from the momentum kick model. The solid, dashed, and dash-dot curves give the total yield, the ridge 
yield, and the pp jet yield, respectively. 

Using the parameters we have extracted from the STAR ridge data as given by (|3ip and ([32|) . we can predict 
the pseudorapidity distribution for the PHOBOS experimental acceptance defined by < 1, < Tytrig < 1.5, and 
0.02 < pt < 2.5 GeV. The theoretical total associated particle yield, which is the sum of the ridge yield and the 
attenuated pp jet yield, is shown as the solid curve in Fig. 5. The theoretical pp jet yield and the ridge component of 
the associated particles are shown as the dash-dot and the dashed curves respectively. The result has been corrected 
for A77 acceptance. The present prediction of the momentum kick model for the near-side jet associated particle yields 
was found to agree well with experimental measurements obtained by the PHOBOS Collaboration [1^ up to large 
IA77I for the region of small pf. 

VII. EXTRACTED INITIAL PARTON MOMENTUM DISTRIBUTION 

It is illuminating to examine the initial parton momentum distribution extracted from the totality of experimental 
data in Figs. 2, 3, and 4. We find that the normalized initial parton momentum distribution at the jet-(medium 
parton) collisions can be represented by dF / dydipptdpt = Aridge(l — x)"" exp{ — -^/m^ + pi) /T} / \/'rn^ ~^ Pt^ where 
X = + p1e^y^~y*' /rrib^ a=0.5, T = 0.5 GeV, and rud = 1 GeV. Here, (y, 0,pt) represent the initial parton 

momentum coordinates. We show explicitly the extracted, normalized initial parton distribution dF / ptdptdy at the 
moment of jet-(mcdium parton) collisions in Fig. 6. It is given as a function of pt for various y in Fig. 6(a), and 
conversely as a function of y for various pt in Fig. 6(b). It has a thermal- like transverse momentum distribution and 
is nearly fiat in rapidity at j/ 0, with sharp kinematic boundaries at large \y\. 




FIG. 6: (Color online) Normalized initial parton momentum distribution dF / dyptdpt extracted from the STAR Collaboration 
data Bail, (a) dF/ dyptdpt as a function of pt for different y, and (b) dF / dyptdpt as a function of y for different pt. 
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In Fig. 6(a), the momentum distribution for y = and high pf has a slope parameter T that is intermediate between 
that of the jet and inclusive bulk particles. This indicates that partons at the moment of jet- (medium parton) collision 
is at an intermediate stage of dynamical equilibration. 

The parton momentum distribution cannot be separated as the product of two independent distributions. The 
momentum distribution as a function of pt depends on the rapidity variable y which affects the boundaries of the 
distribution. The distribution as a function of pt does not change much for y up to y = 2. For y = 3, the maximum 
value of Pt is 1.54 GeV and the distribution changes significantly as the kinematic boundary is approached. For y = 4, 
the boundary of pt is located at 0.55 GcV. 

In Fig. 6(b), the momentum distribution as a function of y for a fixed pt is essentially flat near central rapidities 
and it extends to a maximum value of |y|max that depends on pt, as given by Eq. (|22p . The flat distribution changes 
rather rapidly as it approaches the kinematic limits. The kinematic boundary becomes more restrictive to cover a 
smaller allowed region of y as pt increases. For example, for pt = 0, 1, 2, 3, and 4 GeV, the maximum values of |y| are 
5.36, 3.4, 2.7, 2.33, and 2.05, respectively. The extracted early rapidity distribution exhibits the feature of a plateau 
structure in rapidity. The width of the plateau decreases as pt increases. 

The locations on the kinematic boundaries in Figs. 6(a) and 6(b) depend on the value of yt and mf, which have 
been taken to be yN and niT, respectively in the present analysis. Better determination of these quantities using 
the measured locations of the kinematic boundaries may require more refined measurements of the ridge momentum 
distribution in many locations in pscudorapidity space. 

VIII. EARLY PARTICLE PRODUCTION AT HIGH ENERGIES 

The momentum distribution extracted from the near-side ridge data indicates that the early parton rapidity dis- 
tribution has a plateau structure which extends well into the high rapidity region. The width of the plateau depends 
on Pt. The greater the value of pt, the narrower is the width of the plateau. While the evolution scenario of the early 
rapidity distribution has been outlined in Section VIII of [l^ , we would like to elaborate in more detail the origin of 
the rapidity plateau. 

It should be kept in mind that the plateau rapidity structure has been known in QCD particle production processes 
both experimentally and theoretically. In e~^-e~ annihilation experiments, the produced particles exhibit a rapidity 
plateau structure [H, [13, [H, [s^, H^l . Many earlier theoretical investigations of QCD particle production processes give 
a rapidity plateau distribution when a quark pulls away from an anti-quark at high energies [13, [H, [3, [1^, [1^ . We 
shall review here the theoretical basis for the occurrence of the rapidity plateau in an elementary particle production 
process. 

As an exact solution of particle production at high energies starting from the first principle of QCD is not available, 
many phenomenological models have been presented to describe particle production in nucleus- nucleus collisions ^5^ . 
Common to many of these models (such as the Lund Model, the Dual Parton Model, the Mulitple Collision Model, 
the ART model, the Lexus Model, the Venus Model, and the Glasma Model, ...) is the elementary particle production 
process of a color charge pulling away from an anti-color charge at high energies at the early stage of a nucleus- nucleus 
collision, and the nucleus-nucleus collision consists of many of these elementary production processes. 

We can single out one of the elementary production processes for examination and study the particle production 
process in a model that has many essential features as those in QCD [i^, [4^, [3, [4^, [4y]. The model of QED2 
(42l . [47I . [48l . [49I . [soj are quantum mechanical systems in which a neutral boson exists as a non-perturbativc bound 
state, much as mesons are bound states in QCD. When a positive and negative charge pair are separated in such a 
system, the vacuum is so polarized that the positive and the negative charges are completely screened, in a manner 
similar to the confinement of quarks, in which a quark cannot be isolated. It was demonstrated by Casher, Kogut 
and Susskind p3 | in QED2 that the rapidity distribution of produced boson particles in a system of two oppositely 
charged fermions separating at high relative velocities exhibits a plateau rapidity structure. Such a rapidity plateau 
structure of p roduced particles is indeed observed in high-energy e"*" — e~ annihilation experiments as mentioned 
above (36l . ISTI . [ssI . [sqI . [40 | . The quark fragmentation function obtained from QED2 [53| agrees with that of Field and 
Feynman [58| in their phenomenological treatment of QCD strings. These desirable properties of confinement, charge 
screening, the existence of neutral bound states, and the proper high-energy behavior make it useful for Casher, 
Kogut, Susskind [i^], Bjorken [i^ and many others [H, [i^, [ig to infer the rapidity plateau structure of produced 
particles when a color charge recede away from an anti-color charge at high energies. 

Previously, a scaling argument was presented to reduce QCD at high energies to an effective two-dimensional field 
theory by scaling the longitudinal and temporal coordinates by A and expand the action in power of 1/A [5^. We 
shall try an alternative approach by using the physical argument of transverse confinement to establish the connection 
between QCD and QED2, in order to study particle production in a quantum mechanical framework. 

We consider the elementary particle production process in a flux tube in a nucleus- nucleus collision as a color 
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charge and an anti-color charge separate from each other at high energies. Produced particles are quanta of the 
interacting fields. Depending on the environment temperature, they can be considered as partons in the environment 
of a strongly-coupled quark-gluon plasma and as hadrons in a cold QCD environment at zero temperature. The QCD 
fields inside the tube consists of the gauge fields and the fermion degrees of freedom. At high energies, the gauge fields 
can be greatly simplified by noting that the transverse gauge fields Ax and Ay are expected to be proportional to the 
fermion source transverse velocities, which are smaller as compared to the longitudinal velocity in the z-direction. It 
is reasonable to ignore the transverse gauge fields A^, and Ay so that A^ = {Aq, 0, 0, A^) containing only Ao and A^ 
degrees of freedom confined in the tube. 

The fermion sector can also be approximated. We shall assume that as a result of the non-perturbative non- 
Abelian gauge interaction, transverse confinement is established, and this confinement can be conveniently described 
by a scalar potential m(r) that limits the amplitude of the fermions to the region around the flux tube, as in previous 
descriptions [ssl [60l. [6ll. [6^. [ssj . The Dirac equation for a fermion in the tube in cylindrical coordinates {r,ip,z) is: 

(33) 



{t^'ItTm - eA^) - »7i(r)} *(r, ip, z, t) = 0. 
Following the results of Ref. [gO, [fil] , we seek a solution of the Dirac equation ([33]) in the form 

7''(7r^ - eA^) + m(r) ^p{r,ip, z,t). 



The equation for ijj is 



(34) 



(7^9i + -f'^d2)m{r) 



V'(r,^,z,0 = 0, (35) 



where e is the coupling constant. We note that [a^, Jz] = 0, where Jz = —id/d(j) + <7z/2 is the third component of the 
angular momentum operator. Furthermore, both and commute with the operator acting on ip(r,(l), z,t) in Eq. 
Upon using the representation in Ref. [g^l, the eigenfunction of satisfying a^fj,\ = rixfJ.\ are 



1 I 





1 I 1 



1 



/0\ 
1 I 1 



^0/ V-1/ V / VI/ 

with r]i^2 — +1 and with 773^4 ~ —1. Therefore, we may choose ip^r^ ip, z, t) to be factorized as 



(36) 



'^jAr,f,z,t) ^ ^ fj^^{z,t)Rj^^{r,ip), (37) 

77=-l,l 

with Rj^^ir^ip) to be simultaneous eigenfunctions of Jz and a'^ . The eigenfunctions of Jz satisfying JzRj^r] ~ 
{v + az/2)Rj^^ are 



Rj. i(r, 4>) = 5i.(0e*'^Vi - <?2.(r)e^('^+^) 



M2, 



(38) 



M4- 



(39) 



As a result of the transverse confinement, the gauge fields Aq and Az are confined within the transverse dimensions 
of the flux tube. For high-energy collisions, the transverse dimensions of the flux tube are much smaller than the 
longitudinal dimension. To study the dynamics along the longitudinal direction, it is reasonable to average the gauge 
fields Aq and A^ over the transverse profile of the flux tube. After such a transverse averaging, the dynamics of Aq and 
Az along the longitudinal direction can be approximated to be independent of the transverse coordinates. One can 
then use the method of the separation of variables to separate the equation of motion. By introducing the transverse 
eigenvalue m^, the Dirac equation can be separated into the set of equations in different coordinates. 



[PQ - Ao(z,i)]^ - [pz - eAz{z,t)]^ - m\ - r]ie[dzAoiz,t) - doAz{z,t)]fj^rj{z,t) = 0, 



(40) 



I N .dm{r) 



(41) 
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p\{v + 1) + m^(r) - 



! N .dm{r) 



where 



2 / \ ^ d , d . 

r or or r^ 



(42) 



(43) 



Here, is the eigenvalue for the coupled transverse equations ([4T|) . and (|42p. obtained by imposing the boundary 
condition that the transverse wave functions g\i, and g-2v are transversely confined with a vanishing probability at 
r oo. The eigenvalue mj^ depends on and is independent of the quantum number r]. Some examples of m(r), 
m_L and transverse wave functions have been presented previously (ssl . [6^ . [63j . 

We can write the wave function -0 with the quantum number and a mass mx as a two-component wave function 
in an abstract two-dimensional QED2 space as 



^qod2 



In terms of this wave function, Eq. (|40|) becomes 

7qcd2[Po - e^o(a;\t)] +7qcd2[Pi - e^i(a;\t)] - 
where we re-label the longitudinal z-axis as the x^-axis in QED2, and 



V'qod2(a;\t) = 0, 



^,0 

7qed2 



r 

1 0/ ' 



(44) 



(45) 



(46) 



7qcd2 = *Cr2 = o) ■ 

For brevity of notation, the subscript 'qed2' will be omitted. It should be kept in mind that the transverse state with 
in different transverse excitations correspond to QED2 with different • We shall be interested in the state with 
the lowest mj,. 

The above discussions shows how the fermion and the gauge field in QCD4 in a flux tube can be approximately 
mapped into elements in QED2 for high-energy processes. Although the non-Abelian nature of the gauge field in 
QCD is needed to lead to the formation of the confining flux tube, the non-Abelian property is not needed for particle 
production in QED2 at high-energies. An Abclian QED2 fleld theory possesses the desirable properties of confinement 
and charge screening, and it suffices to describe the particle production process at high energies. Furthermore, in the 
non-Abelian field tensor 

Fl^ = aoAi - ai Af, + gr^'KAX. (48) 

the non-linear quadratic term contains the product of Aq and A^. One can convenient choose the Coulomb gauge 
= such that the non-linear quadratic term does not contribute. We can therefore ignore the non-Abelian nature 
of the gauge fields and approximate them to be 

= dt^A, - d,A^, (49) 
where ^,,v = The fermions give rise to the current 

f = e^pjf'i), (50) 

which generates the gauge fields according to 

d^F'"' = -j". (51) 

Eqs. (|45)) . (|49|) . and ((5T|) constitute the equations for the quantum mechanical system of QED2 with a fermion of 
mass m±. Thus, by assuming QCD confined within a flux tube, the longitudinal dynamics of the system can be 
approximated as those of QED2 with a mass m±. The gauge fields (/x = 0, 1) depend on the fermion field tp. The 
fermion field -0, in turn, depends on the gauge field A^. The coupling is quite complicated and leads to a non-linear 
problem of great complexity. Remarkably, Schwinger found that QED2 involving massless fermions with the gauge 
interaction is equivalent to a free boson field with a mass /xq = e/y^, where e is the coupling constant [i?} . 
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IX. PARTICLE PRODUCTION AS AN INITIAL- VALUE PROBLEM IN BOSONIZED QED2 

In mapping elements of QCD4 approximately into elements of massive QED2, what is the relationship between 
the coupling constant g in QCD4 and the coupling constant e in QED2? By limiting the motion and the source 
distribution to reside in the longitudinal direction, the coupling constant e in QED2 acquire the dimension of a mass. 
The confinement property is a non-perturbative property of QCD4. The coupling constant e in QED2 should therefore 
be non-perturbatively related to g. The relationship can be retrieved by comparing non-perturbative quantities. In 
QED2 with the Coulomb gauge Ai — 0, the interaction energy between a quark and an anti-quark separated at a 
separation of x is e^Ao(a;) = e^x/2. On the other hand, in QCD4, the non-perturbative confining interaction energy 
between a quark and an anti-quark is bx where b is the string tension. Therefore, equating the two interaction energies, 
we find a relation between e in QED2 and the non-perturbative string tension b in QCD4, 

e = V2b. (52) 

If we take the string tension to be 6 = 1 GcV/fm, then e = 0.628 GeV. The boson mass in massless QED2 is 
^0 e/V^ 0.354 GeV. 

The case of massive QED2 with a fermion mass can be studied by bonsonization. It is equivalent to the system 
of free bosons of mass interacting with an interaction that depends on [iol . [50| . Does QGD corresponds to 
the case of strong coupling with e 3> or the case weak coupling with e <C m±l The case of strong coupling is 
characterized by a quasi-free bosons with confining fermions and charge screening, while the limit of weak coupling 
approaches free Dirac theory with almost free fermions dressed up as bosons having a mass close to the free fermion 
rest mass (50| . 

We can estimate m±_ to be of the order of ?i/(tube radius) where the flux tube radius is of order 1 fm, leading to 
m_L ~ 0.2 GeV. We have /xq S> ?7i_l which correspond to the case of strong coupling with fermion confinement and 
color-charge screening, rather than quasi-free Dirac particles. Accordingly, the mass-perturbation theory can be used 
to discuss the particle production process in our case of massive QED2. 

In the mass-perturbation theory, the unperturbed theory is massless QED2 and the mass is treated as a 
perturbation. Up to the second order in m^, the mass perturbation theory gives a quasi-free boson with a mass M 
given by [i^ 

M2 = + 2e^^iQmx_ + l.OeTSe^Tmi, (53) 

where 7 = 0.5772 is the Euler constant. We therefore have M ~ /zo + e^m±. For our case of /xq = 0.354 GeV and 
m_L ~ 0.2 GeV, we get M ~ 0.71 GeV, which comes close to the spin-spin averaged mass of 0.62 GeV for the -K-p 
pair. Thus, the boson of massive QED2 finds its correspondence as the boson in QGD that splits into tt and p when 
the spin-spin interaction is taken into account. 

As mass perturbation theory is based on massless QED2 with toj_ as a perturbation, the application of the theory 
to particle production process involves in using the results of massless QED2 and replacing the boson mass pq in 
these massless QED2 results by the corrected mass M . In practical applications, this amounts to replacing p^ with 
the physical mass, including the effects of the effective mass increase due to the transverse momentum. As pions 
are the most predominantly produced particle, the phenomcnological treatment then involves in replacing p^ by 

We can review the rapidity distribution for massless QEp2 obtained previously [1^. The relation between the 
bosonic and the fermionic quantities in massless QED2 is ji^ . |49| . [soj 

f = -ee^^d^cj^/^, (54) 

where is the fermionic current which can be taken to be a real quantity, and e^" is the antisymmetric tensor 
= — £01 = — 1- We note that, as is a vector field and e^^ is a pseudotensor, the field is a real pseudoscalar 
field, and it represents the color electric field F"^, as F^^ = e0/y^. If the current arising from the fermions is 
initially known, then the dynamics of the pseudoscalar field 4> can be inferred at all times. Treating the problem as a 
system of quasi-free bosons with a mass M, the initial value conditions will allow us to determine the dynamics of the 
system. To apply the results to our case, we will work within mass perturbation theory which is a quasi-free boson 
system with e/y^ replaced by M . 

Given an initial fermion charge distribution j^(^x, t = 0), its Fourier transform is i^{p^) is 

.7>(/) = ^ 1 dxe-^P'^f{x, 0). (55) 



16 



We show previously [44| that the momentum distribution of the bosons is then given by 



dN 



dpi 2p0e2 

and the rapidity distribution of the produced particles is 

dN 
dy 



P i IP 



TT 

2? 



^f{p')+jHp') 



P 



^f{-p')+j'{-p') 



(56) 



(57) 



This gives a simple relation between the rapidity distribution and the Fourier transforms of the initial fermionic charge 
current. 

We can review how this initial-value problem in massless QED2 can be formulated for the case of a positive charge 
ve separating from a negative charge — i^e with a center-of-mass energy [31 ■ We work in the center-of-mass system 
and start at t = with the charge and anti-charge pair superimposed so that the total charge density of the system 
at < = is zero: 



/(.T,0)=0. 

To construct the initial longitudinal current, we introduce a distribution that depend on a. 



j\x,t) = - {^[tanh((x + t)/<j) + 1] + ^[tanh((x - t)/a) + 1]} 



(58) 



(59) 



In this case, the initial current which arises from a charge ve moving in the positive x direction and another charge 
—lye moving in the negative x direction is given by 



(7 cosh (x/a) 



(60) 



In the limit as a approaches zero, the above current is proportional to a delta function. The diffusivity a is related 
to the total invariant mass ^/s of the system; using the energy P'^ = y/s at the initial time t = 0, we obtain a relation 
between a and \/s: 



For this current distribution j^{x, 0) the Fourier transform of j^{x, 0) is 



3'{p') = - 

and the rapidity distribution in massless QED2 is [44| 

dN 



/27r sinh(7rpiCT/2) ' 



dy sinh^ ^ ' 



where 



i^^TT^/iQ sinhy 



(61) 



(62) 



(63) 



(64) 



The rapidity distribution therefore shows a plateau structure around y ~ 0. In the limit of very high energy, the 
rapidity distribution is dN/dy = , which agrees with the result of Casher [i^ (for v = \). 

Within the mass perturbation theory, we can approximate the particle production process of massive QED2 using 
the results from massless QED2 and replacing the /ig of massless QED2 in Eq. by M in massive QED2, with the 
result 



V -K M sinhy 



(65) 



Using Eqs. ([63]) and ([65]) by replacing //q with the phenomenological mass M = \J'm\ -I- (p_L,7r)^ = 0.30 GeV and 
V = 2.45 gives a good phenomenological fit to the d7Vjr± /dy data in e+-e~ experiment at -^s = 29 GeV (in Fig. (40) 
of P). 
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X. EVOLUTION OF THE MEDIUM PARTON MOMENTUM DISTRIBUTION 

We conclude from our discussions in the last two sections that in addition to experimental evidences for the rapidity 
plateau in elementary QCD particle production processes, theoretical investigations in plausible models also show the 
occurrence of a rapidity plateau when a color charge pulls away from an anti-color charge at high energies. As a 
nucleus-nucleus consists of elementary production processes of string fragmentation, there can be a similar plateau 
structure in the rapidity distribution of the produced medium partons, consistent with the parton rapidity plateau 
we have extracted at the early stage of the nucleus-nucleus collision. 

In nucleus- nucleus collisions, this early parton momentum distribution can be probed by a jet produced in the early 
stage of the collision. Those medium partons kicked by the jet subsequently materialize as ridge particles and they 
retain the property of the rapidity plateau. 

The plateau rapidity structure of the early parton momentum distribution differ from the Gaussian rapidity distri- 
bution of the bulk matter [g^, [H, [63| • How does one understand such a difference? 

It is important to point out that the early momentum distribution represents the momentum distribution at the 
early stage of the nucleus-nucleus collision as it involves the direct reaction with the jet, which occurs only at the early 
stage of the nucleus-nucleus collision. On the other hand, the momentum distribution of the bulk matter represent 
the momentum distribution of the bulk matter at the end-point of the nucleus-nucleus collision. A considerable period 
of time separates the beginning, early stage of the nucleus-nucleus collision and the end-point of the nucleus-nucleus 
collision. Significant dynamical evolution must have occurred between these two separate time points, as described 
schematically in Fig. 9 of Ref. [l9| . The time evolution of the momentum distribution will make the end-point 
momentum distribution of the bulk matter different from the early parton momentum distribution. 

Evidence for the occurrence of a dynamical evolution of the momentum distribution presents itself in the difference 
of (i) the transverse momentum distribution extracted at the moment of the jct-(mcdium parton) collisions, and (ii) the 
transverse momentum distribution of the bulk matter at the end-point of the nucleus-nucleus collision. The extracted 
early parton transverse momentum distribution, as given by Eqs. (jl9p . and (|32p . has a thermal-like distributions, 
with an initial inverse slope T = 0.5 GeV that is slightly greater than the inverse slope of the end-point transverse 
momentum distribution, consistent with the direction of transverse momentum evolution from a higher inverse slope 
T to a lower inverse slope T value Q. We expect that the rapidity distribution will likewise evolve and its shape 
will change with time. There is no reason to expect that the longitudinal momentum distribution at the early stage 
of the nucleus-nucleus collision should be the same as the corresponding longitudinal momentum distribution at the 
end-point of the nucleus-nucleus collision. 

To understand the evolution of the medium parton momentum distributions, we should think of the full momentum 
distribution to be a six-dimension distribution function F{t, p, t) of the medium that depends both on the spatial 
and momentum coordinates, as well as on the lapsed time t. The parton momentum distribution extracted here is 
in effect an average of this the six-dimensional distribution function F over spatial and temporal coordinates of the 
collision points in the early stage of the nucleus-nucleus collision, using the jet as a probe. After the early stage 
of jet-(mcdium parton) collisions, partons from one position will collide with partons of adjacent positions. These 
collisions will modify the momenta of the colliding partons, leading to a change of the distribution function F{v, p, t) 
as a function of time. How the evolution will take place is a problem of great complexity that depends on models with 
many unknown theoretical elements [s^. Nevertheless, one expects that starting with a non-isotropic plateau rapidity 
distribution that is much elongated in the longitudinal direction, a collision of two partons with large and opposing 
longitudinal momenta in adjacent spatial locations will redistribute the partons from the longitudinal direction towards 
the transverse directions, with a decrease in the longitudinal momenta of the colliding partons. Hence, the evolution 
will smooth out the anisotropic plateau rapidity structure to a significant degree as time proceeds. 

XI. DEPENDENCE OF THE FRAGMENTATION FUNCTION ON THE JET- (MEDIUM PARTON) 

COLLISION NUMBER 

We turn now to investigate the geometry-dependent part of the ridge and trigger particle yields, as given previously 
in Eqs. ([7]), and (fT5|) . The jet fragmentation function in these equations depends on the number of collisions 

N (or Nk) suffered by the jet parton, and the observed physical quantities depend on the distribution P{N). We 
envisage jet partons to be produced by binary nucleon-nucleon hard-scattering processes and we focus our attention 
on one of the jet partons. We consider the jet parton to occur at ryjet = such that the jet momentum pj, the trigger 
particle momentum ptrig, and the momentum kick q all lie in the mid-rapidity transverse plane in the collider system, 
pointing in the same direction. The vector symbol for these quantities can be understood. 

We envisage that in the passage of the parent jet parton in the dense medium, the jet parton with initial momentum 
Pj imparts a momentum q to each kicked medium parton, and loses a momentum Nq after N jet-(medium parton) 
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collisions. In addition to coUisional momentum loss, the jet parton can lose momentum by gluon radiation As 
the radiated gluon will likely come out in a cone along the jet direction in random azimuthal angles, the average 
momentum loss due to gluon radiation (A,,) lies along the jet direction ej. We can parametrize the radiative gluon 
momentum loss phenomenologically by |Ar| = Nqr where the qr value obtained from experimental data will need 
to be compared with theoretical models. In practice, the coUisional and radiative momentum losses appear together 
as the sum total (q + qr) in the fragmentation frmction [Eq. (j7ip ]. Furthermore, there can be additional attenuation 
Ca due to absorptive inelastic processes of removing the jet from the jet channel. Only the sum of the absorptive, 
coUisional, and radiative contributions, leading to the total attenuation coefficient can be obtained by comparison 
with experimental jet quenching data [Eqs. (I73|) and ([75]) ]. 

Upon including the momentum loss due to jet- (medium parton) collisions and gluon radiation in the momentum 
kick model, Eq. ([7]) becomes 

f dN ^"""^ 

iVtHg= / dp,—^ P{N)e-^'^''D(j>,,,^-p,-N{q + qr)). (66) 

We wish to write out the dependence of the fragmentation function Z? on iV in the above equation explicitly. The 
dominant contribution of the jet production process comes from gluon-gluon collisions [68j . The relevant fragmentation 
function of fragmenting a pion out of a gluon at the momentum scale Qq, can be written in the form [68| 

zD[z,Ql)^C^{l-z)'^\ (67) 

where z ~ Pti-ig/Pj and is a constant. In perturbative QCD, the fragmentation function near z = I varies with the 
QCD momentum scale Q according to [H, [6^ 

zD{z,Q') zD{z, Ql)e'-''^\- In . (68) 

where G = 4/25, and s = ln[ln((5^/A^)/ ln((5Q/A^)]. After the jet suffers N jet-(mcdium parton) collisions with 
the medium partons, the fragmentation function for a pion to fragment out of the final jet of momentum zn = 
Ptrig/[Pj - N{q + qr)] is 

D{zN, 0') = ae°-«9^^"exp{- Inz^ + ai ln(l - zjv) + 4Gsln[(- InzAr)]}— (69) 

I (ai + 1 + AGs) 

Upon expanding the exponent index of the above function in power of N{q + qr)/Pj and retaining the term first order 
in N{q + qr)/pj, we obtain the dependence of the jet fragmentation function on the jet-(medium parton) collision 
number N , 

D[ptriJ{Pj - N{q + qr)),Q^) ~ I?(ptrigM-,Q')e^«"^, (70) 

where 



Ol ^ AGs Pj \ Ptrig(g + Qr) 

Ptrig/Pj ' 1-ptrig/Pi ln(pj/ptrig)ptrig/ 



^D=[z — + -, 1 — + TZ7Z-7Z 12 ■ ("1) 



Substituting this relationship into Eq. ([66|) . we get 

^tog= / dp,— li?(ptHg; P,) J2 (72) 

P-' N=0 



where we have combined (a with (^u, 



C-Ca+Cc (73) 



Because of the normalization condition Eq. ([5]) for the fragmentation function and the definition of dNj/dpj in Eq. 
(H]), we obtain 



N„ 



Ntris = Nun J2 PiN)e~^'', (74) 



Af=0 
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and we get the jet quenching measure 

R^^ = pE^Y.P^N)e~CN^ (75) 

JVbin 

and the average number of jet- (medium parton) coUisions per trigger 

(iV) = (Nk) = J2 NP{N)e-^^ / J2 PWe-'^^. (76) 

N=l ' N=0 

The presence of the attenuation factor e"*"^ impUes that detected trigger particles are hkely to originate near 
the surface where the number of jet- (medium parton) collisions N is smallest. The quantity (^a is not known. The 
quantity depends on g -I- and pj . We can estimate the contribution of coUisional contribution to the value of (^d 
by inferring the approximate average value of pj for ptrig ''^ 5 GeV. As pj ^ (ptrig + {N)q + 3iVjot?jot) with (N) ~ 6, 
q ^ 1 GeV, iVjet=0.75, Tjet = 0.55 GeV, we estimate that pj - 2.5ptrig- If we use = p], Ql = 3 GeV^ ai = 1.5 
and A = 0.5 GeV as in [6^ we can use Eq. ((7T|) and estimate the contribution to from q to be approximately 0.22. 
There can be additional contributions from the radiative energy loss qr- We shall set C as a free parameter to describe 
the experimental Raa data by searching for ( around the neighborhood of about 0.22. We find in Section XIII that 
the experimental jet quenching and ridge yield data are consistent with a value of C = 0.20 which comes very close to 
the value of 0.22 estimated here. 



XII. GEOMETRY DEPENDENCE OF TRIGGER AND RIDGE YIELDS 



Because the ridge particle yield has been measured on the basis of the yield per trigger particle, it is necessary to 
determine the trigger yield A'trig as a function of the geometrical variables, in addition to determining the number of 
ridge particles. The trigger particle yield is quenched due to the energy loss of the jet parton. We therefore need to 
study jet quenching and follow the trajectory of the jet. 

From our earlier considerations, the relevant physical quantities are given in terms of N, the number of medium 
partons kicked by the jet on its way to emerge from the medium. This quantity N, in turn, is equal to the number of 
jet-(medium parton) collisions Nk suffered by the jet parton. We assume for simplicity that the energetic jet parton 
travels along a straight line trajectory with a velocity nearly the speed of light, making an angle (j)s with respect to 
the reaction plane, 0^ = (/)jot — 0Rp. Using the mid-point O between the centers of the two nuclei as the origin, we 
set up a transverse coordinate system for the jet source point bo and the jet trajectory point b' as shown in Fig. 7. 




FIG. 7: (Color online) The transverse coordinate system used for the jet source point bo and the trajectory point b'. The 
coordinate origin is located at the midpoint O between the two colliding nuclei, whose centers are located at Oa and Ob 
separated by an impact parameter b . The jet trajectory lies along I and makes an angle (ps with respect to the reaction plane. 

We consider the jet source point at bo, from which a mid-rapidity jet parton originates. The number of jet- (medium 
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parton) collisions along the jet trajectory making an angle (f>g with respect to the reaction plane is 

iV,(bo, 0.) = J\dl^^^^ [b'(bo, 0.)] , (77) 

where < ^ < oo parametrizes the jet trajectory, a is the jet-(medium parton) scattering cross section, and 
dNpaTton{b')/dV is the parton density of the medium at b' along the trajectory Jet- (medium parton) collisions take 
place along different parts of the trajectory at different I and involve the medium at different stages of the expansion. 
They depend on the space-time dynamics of the jet and the medium. To follow the jet- (medium parton) collisions 
along the jet trajectory, we need a time clock to track the coordinates of the jet and the motion of the medium. We 
start the time clock for time measurement at the moment of maximum overlap of the colliding nuclei, and the jet is 
produced by nucleon-nucleon collisions at a time t ^ fi,/(10GeV) which can be taken to be ^ 0. The trajectory path 
length / is then a measure of the time coordinate, t ^ I, which is needed to follow the longitudinal and transverse 
expansions of the medium. 

The trajectory point b' depends on the jet source origin point bo and the jet azimuthal angle 0s as 

b'(bo,0s) = (6^,&p = {bax + lcos(j)s,boy + lsm(j),). (78) 

If we approximately represent the modification of the density arsing from longitudinal and transverse expansion by 
an effective time parameter using the initial parton density in the following approximation, 

N,iho,cj^s) = — r dl'^^^^ [b'(bo,0.)]. (79) 

then we find that the data can be described by C = 0.22 and a Jt^s ^ 0.025 fm. Such a picture only gives a rather 
crude description of the path-length dependence of the ridge yield. 

In order to give a more realistic picture, we describe the medium by an expanding fluid with an initial density given 
by the distribution of the participants at the moment of maximum nuclear overlap. We assume that the longitudinal 
expansion begins at the moment of maximum overlap as the initial momenta are directed along the longitudinal 
direction. A period of time to is however needed to convert the longitudinal kinetic energy into entropy to produce 
particles with a transverse mass. The transverse hydrodynamical expansion can then commence at t > ^o- The time 
for producing a particle with a typical transverse mass of about 0.35 GeV is fi./(0.35 GeV) ^ 0.6 fm/c, which is also 
the time estimated for the thermalization of the produced matter [t^I- We therefore take tg = 0.6 fm/c. 

As we will focus our attention in the mid-rapidity region where experimental data are available, Bjorken hydrody- 
namics and Landau hydrodynamics (tT], [T^I coincide [zl, [tI] and we can use Bjorken hydrodynamics to describe 
the longitudinal expansion. For a hydrodynamical system undergoing Bjorken longitudinal expansion, the transverse 
expansion can be described by the hydrodynamical solution of Baym et al. (75| . Using the method of characteristics, 
they find that the energy density and velocity field in the transverse direction can be described well approximately 
by analytic formulas. 

Accordingly, we follow the jet along its trajectory at I specified by the trajectory point b' with a transverse magnitude 
b' = |b'| (measured from the origin O) at at time t = I. The time after the onset of the transverse expansion is then 
tf( = t — to, and a rarefaction wave travels from the transverse radius R inward with the speed of sound Cg- The 
dynamics is different whether the rarefaction wave has reached this medium point b' or not. The transverse space of 
the medium can be divided into Region I and II. 

In Region I characterized by &' < i? — Cgtji, the rarefaction wave has not reached this medium point at b'. In this 
region, the medium has not started to expand transversely with transverse velocity v± = while the longitudinal 
expansion has already commenced. Due to the longitudinal expansion the density is depleted and the temperature is 
decreased as [zl] 

Tcx(io/t)"'. (80) 
As the entropy density and number and entropy densities are proportional to T^/^^ , we have 

-^y- ' ^) = (^init ,t = to)j, (81) 

with 6Jj^;j = b' in Region I. 

In Region II defined by R — Cstfj < b' < R + tn, where the inward-traveling rarefaction wave has passed through 
already. The medium is expanding transversely outward and the transverse velocity v± at the point b' at the time t 
is 

b' -R+ CstR 



tR + Cs{b'~R)' 



(82) 
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The transverse velocity v± is unity (speed of light) at the surface point R + tji, and is zero at at the point b' — c^tfi 
where the rarefaction wave has just arrived. The medium temperature in this region is given by [75| 

m<).w;,„.,<.,.,(|^|^)■■"(^f)^ (83) 

where h[^^^ = b' — v±tR is the initial position that reaches b' at tji, and if b' — v±tii < we set b[^-^^ = 0. Here the 
exponential index d is |75j 



c2 



1 

1 



l-v^ib',t)c. 



(84) 



The corresponding medium number density in the coordinate frame in which the trajectory I is measured (Fig. 7) is 
therefore 

(b,t)=j (&i„it, t = to) ( -^—T, E T—] [-] ' (85) 

where 



dV ' ' dV ' "'\tji + b'-R l + csj \t 



7=-=^, (86) 

and 7 is to take into account the change in the medium number density due to the flow velocity of the medium along 
the transverse direction. While the boundary R is independent of the azimuthal angle of b' for the central collision, 
the boundary radius is a function of the azimuthal angle for non-central collisions. Wc shall assume that the relations 
between the density and the radius given above remain applicable by using a radius R that depends on the azimuthal 
angle. These results of the number density at various transverse points allow one to obtain the absorption exponent 
index for a jet to pass through an expanding medium. In numerical calculations, we take the speed of sound to be 
cs = 1/V3. 

The medium parton density dNpmton/dV at {b[^i^,t = to) is related to the parton transverse density dA^parton/c^b 
at to by 

^^^ibL.t = to) - ^^^ibUt = to) (87) 

We can relate the initial parton number transverse density dA'^parton/c'b' at t = to to the corresponding participant 
initial number transverse density dA'part/rfb' as 

C^-^parton _ C^-^parton C^-^part _ C?-^part 
db' ~ diVpart dh' ~ dh' ' 

where K=dA'parton/diVpart is the number of partons per participant. A previous collection of data gives 
iVch/(iVpart/2) = 28 for ^/sWn = 200 GeV and 16 for ^/FWn = 62 GeV (see [zi, ^). If we use the parton-hadron 
duality and count the parton number by (3/2) times the charged multiplicity of detected hadrons, then we get 

^ ^ f21 for ^i]^ = 200 GeV 

[12 for = 62 GcV. ^ ' 

A given source point bo and a given azimuthal angle c/js will lead to iVfc(bo, (fis) number of kicked medium partons, 
which we shall identify by parton-hadron duality as ridge particles. The jet number transverse density is given by the 
binary nucleon-nucleon collision number transverse density, as nucleon-nucleon collisions are the source of jets. We 
need to weight the number of kicked medium particles by the local binary collision number element dbo x dNhin/dho. 
The normalized probability distribution P(iV, 0^) with respect to the number of ridge particles (or jet-(medium parton) 
collisions) is 

P(iV,0,) = ^ / dbo^(bo).5[iV-7Vfc(bo,(A.)], (90) 

^Vbin J "Do 

which leads to the desired normalization of the distribution P{N, (ps), 

dNPiN,^s)^h (91) 
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and the total number of binary nucleon-nucleon collisions 

iVM„=/dbo^. (92) 

J "Do 

Thus, the number of ridge particle yield per trigger particle (or the number of jet- (medium parton) collisions per 
trigger) at an azimuthal angle (ps, averaged over all source points of binary collisions at all bo points, is 

iVfe(0.) = / NP{N, (t>s)e-^'' dN / [ P{N, 4>s)e-^'' dN. (93) 



In practical calculations, it is convenient to discretize N by replacing the delta function in Eq. ((90)l as 

6{N ~ Nk) ^ {e[Nk ~{N~ AN/2)] - e[Nk - {N + AN/2)]} /AN. (94) 
Upon choosing AA^ = 1, the normalization condition ([91]) becomes 

£p(iV,0,) = l, (95) 

as in our previous definition with the 0s dependence now explicitly written out. Equation (j90p in the discretized form 
of N becomes 

P(7V,0,) = -^ / dbo^^(bo){e[(iV-A7V/2)-iV,.(bo,0,)]-e[(iV + AiV/2)-iVfe(bo,0,)]}/A7V. (96) 
A^bin J dho 

This equation facilitates the evaluation of P{N, (pg). For a given cpg, we evaluate Nk{ho, <f)s) at different source points 
bo, place the quantity dNhin/dho at the appropriate [N — AN/2 < Nk{ho,<j>s) < N + AN/2] bin, and accumulate 
the contributions from all jet source points at all bo. The accumulated distribution, divided by A^binAA^ is then the 
distribution function P{N, (ps). For these calculations, we need the transverse densities of the binary nucleon-nucleon 
collisions and participant numbers. The transverse density of binary nucleon-nucleon collisions in Eq. (|96p can be 
obtained from the Glauber model to be 

^(bo) = ABr(b^o)r(bBo)^ir , (97) 
dbo 

where b^o = bo + b/2 and b^o = bo — b/2. The quantity a^-J^ is the nucleon-nucleon inelastic cross section, which 
we can take to be 42 mb at ^/sJ/n — 200 [tHI- The participant number transverse density needed in Eqs. ([79]) and 
([88]) along the jet trajectory can be similarly obtained from the Glauber model to be 

^^ih')=.AT{hU) + BTih's), (98) 

where the transverse coordinates are given by b^ = b' + b/2, and h'g = b' — b/2. These relations allows us to use Eq. 
([96)1 to evaluate P(Af, (/)s). The distribution P{N,(t)s) can then be used in Eq. (|75p to evaluate Raa{<Ps), and in Eq. 
([75)) or ([M)) to evaluate Nk{4>s)- After Nk{4>s) and Raa{4>s) have been evaluated, we can average over all azimuthal 
angles (pg and we obtain the ridge particles [or jet-(medium parton) collision] per trigger 

(A^,) = / dcPsNk{cj)g)/{TT/2), (99) 
Jo 

and 

{Raa) = / dcPsRAAicps) /{7t/2), (100) 
Jo 

which is usually expressed just as Raa- 

Eq. ([15|) separates the ridge particle distribution into a geometry-dependent part, {fR){2/3){Nk), and the normalized 
ridge momentum distribution. dF/dp. From the magnitude of the ridge yield, we have extracted phenomenologically 
in Section VI the values of {f]i) (Nk) — 3.8 for central AuAu collisions at ^/snn = 200 GeV. Ridge particles after 
production are attenuated before reaching the detector. It is reasonable to take the average ridge particle attenuation 
factor (/fl) to be the same as the average attenuation factor for jet component particles, fj = 0.632, as both types 
of particles come out from the interacting region. We then get an estimate of (A^) = {Nk) = 6.0 as the total number 
of kicked partons per trigger for the most-central AuAu collisions at 200 GeV. For numerical purposes, we shall use 
these average numbers as references, keeping in mind however that they depend on the attenuation factor {ffi} that 
may be uncertain. 
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FIG. 8: (Color online) The ratio Raa of the high-pt tt" trigger yield in AA collision, as a function of the number of participants 
A^part- The solid curve is the theoretical result from the momentum kick model, using ^ = 0.20 and a — 0.7 mb. The data 
points are from PHENIX high-pt tto measurements for Au Au collisions at y^sjvjv = 200 GeV [t^ . 



XIII. COMPARISON OF RIDGE YIELD WITH EXPERIMENTAL CENTRALITY DEPENDENCE 



For a given impact parameter and azimuthal angle (j)s, the unknown parameters are C, and a. Although all quantities 
depend on these two parameters, the quantity Raa for the quenching of the trigger is more sensitive to C, and the ridge 
yield per trigger is more sensitive to a. We find that the totality of experimental data of the ccntrality dependence of 
Raa and the centrality dependence of the ridge yield, can be explained well when we use 

C 0.20, and <t = 1.4 mb. (101) 




FIG. 9: (Color online) The ridge yield per trigger as a function of the participant number A'part for nucleus-nucleus collisions 
at -^sjvAr=200 GeV. The solid curve gives the theoretical result for AuAu collisions in the momentum kick model. The solid 
circular data points are from the STAR Collaboration 0] . 



We discuss here the comparison of theoretical results with experimental data using the above two parameters for 
AuAu collisions at y's^v Ar=200 GeV. Solid circular points in Fig. 8 give experimental PHENIX Raa data for high-pt 
TT^ yields [t^. Theoretical Raa result in the momentum kick model obtained with Eq. (|75p as a function of the 
participant number is shown as the solid curve. It gives good agreement with experimental Raa data. The quenching 
of the jet is well accounted for in the momentum kick model. 
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The STAR AuAu ridge yield per trigger at y^sjvw=200 GcV, shown as sohd circular points in Fig. 9, are taken 
from Fig. 2 of 0. They were obtained for 3 < pt,trig < 4 GeV and 2.0 < pt, associated < Pt.trig- The sohd curve in 
Fig. 9 is the theoretical ridge yield per trigger for AuAu collisions. It has been normalized to match the data point 
(within errors) for the most-central collision examined in Sections IV- VI. Our comparison of momentum kick model 
results and the experimental data in Fig. 9 indicates that theoretical ridge yields per trigger agree with experiment. 
It increases as the number of participants increases. 

The value of C = 0.20 is nearly the same as our earlier estimate of C = 0.22 arising from coUisional jet momentum loss 
alone. This indicates that coUisional momentum loss may contribute the dominant component of the jet momentum 
loss, but more research on theoretical predictions for arc needed to separate out the different absorptive and radiative 
contributions. The cross section corresponds to a parton interacting radius of 0.21 fm. which means that a parton 
having the entropy content of a hadron appears to the jet probe as a strongly interacting scattering disk with a radius 
of 0.21 fm. 



XIV. DEPENDENCE OF THE RIDGE YIELD ON COLLIDING NUCLEI MASSES AND ENERGIES 



Whereas our attention so far has been focused on AuAu collisions at ^/nn = 200 GeV, we would like to investigate 
in this section how the ridge yield scales with the mass numbers and the energies of the colliding nuclei. Experimental 
data for such an analysis have been obtained by the STAR Collaboration with the acceptance region of 3 < _pt,trig < 4 
GeV and 2.0 < pt. associated < Pt.tiig [1, This region of acceptance is slightly different from the acceptance region 
for Fig. 9 used in Ref. The measurements of the ridge yield for AuAu and CuCu collisions at ^/snn =200 and 62 
GeV within the same acceptance region in [1, allows a consistent comparison across mass numbers and energies 
of the colliding nuclei. The experimental ridge yield as a function of the participant numbers are shown in Fig. 10 
as solid points for AuAu collisions and open points for CuCu collisions [1, The circular data points arc for 

i/sjvjv = 200 GeV and the square points are for 62 GeV. One notes that the ridge yield appears to increase with 
increasing number of participants and increasing collision energies. The ridge yield for CuCu collisions is small and 
contains large systematic errors. 
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• STAR AuAu 200 GeV data (Bielcikova et al.) 

O STAR CuCu 200 GeV data (Nattrass et al.) 

■ STAR AuAu 62 GeV data (Nattrass et al.) 

□ STAR CuCu 62 GeV data (Nattrass et al.) 
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FIG. 10: (Color online) The ridge yield per trigger as a function of the participant number A'part for nucleus-nucleus collisions 
at y^ijvjv=200 and 62 GeV. The solid curves are theoretical results for AuAu collisions and the dashed curves for CuCu collisions 
in the momentum kick model. The solid points represent AuAu data and the open points represent CuCu data, from the STAR 
Collaboration 



We show in Fig. 10 the theoretical ridge yield for AuAu and CuCu collisions as a function of the number of 
participants for ^ys^N = 200 and 62 GcV. The experimental and theoretical ridge yields are matched for the most- 
central AuAu collision data point at y^sjvw =200 GcV. Our comparison of momentum kick model results and the 
experimental data at different energies, different nuclear masses, and different participant numbers indicates that the 
theoretical ridge yield agrees well with experiment. For the same nucleus-nucleus collision at different energies, the 
theoretical ridge yield scales approximately with k, the number of medium partons produced per participant, which 
increases with increasing collision energy as (In-^s)^ (76j . For the same collision energy, the theoretical ridge yields 
per trigger for CuCu collisions follow approximately those of the ridge yields for AuAu collisions, when plotted in 
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terms of the number of participants. 

XV. DISCUSSIONS AND CONCLUSIONS 

The experimental near-side ridge data have guided us to the momentum kick model as a description of the ridge 
phenomenon. The narrow cone of associated particles along the jet direction reveals that the trigger particle is 
connected with the occurrence of a jet. The yield of the associated particles as increasing with increasing participant 
numbers and the similarity of their inverse slope reveal that the ridge particles come from medium particles. The 
short-range behavior of the strong interaction due to color screening and the narrow azimuthal correlation of a ridge 
particle with the trigger particle reveal further that the ridge particles and the jet are related by collisions. Hence, a 
picture of the momentum kick model emerges as a plausible description of the ridge phenomenon. 

In the momentum kick model, a jet parton produced in high-energy heavy-ion collisions makes collisions with 
medium partons. The kicked medium partons subsequently materialize as ridge particles while the jet loses energies 
and fragments into the trigger particle and other fragmentation products. 

The implementation of the momentum kick model can proceed numerically by a Monte Carlo approach, following 
the trajectory of the jet and the medium particles as the medium evolves. The space-time dynamics of the medium 
and the jet is a problem of great complexity and contains many complex, unknown, and non-perturbative elements. 
However, before we implement such an elaborate undertaking, it is useful to explore with simplifying assumptions 
whether the momentum kick model contains promising degrees of freedom. 

Following the dynamics of a jet and and its interaction with the medium, we show how the ridge yield can be 
greatly simplified by using the average distribution of the medium particles and the average momentum kick. We are 
then able to separate the ridge particle yield into a factor which depends on the average number of partons kicked 
by the jet and another factor related to the (average) momentum distribution of the kicked parton after acquiring a 
momentum kick from the jet. The ridge particles therefore carries information on the momentum distribution of the 
partons at the moment of jct-(medium parton) collisions. They also carry information on the (average) magnitude 
of the momentum kick a medium parton acquires. These complications of space-time dynamics of medium and jet 
partons have been subsumed under the probability distribution P/v(iV), which depends on geometry, medium parton 
dynamics, jet parton trajectories, and jet- (medium parton) cross sections. 

The medium partons kicked by the jet materialize as ridge particles can be used to extract the early parton 
momentum distribution. The extracted early parton momentum distribution provides valuable information for the 
mechanism of early parton production and the later evolution of the system toward the state of quark-gluon plasma. 
For central AuAu collision at y^s^vw — 200 GcV, we find the extracted early parton momentum distribution to 
have a thermal-like transverse distribution but a rapidity plateau structure whose width decreases as the transverse 
momentum increases. We should note that plateau rapidity structure has been known in QCD particle production 
experiments [H, [13, [H, [H, and in QCD particle production theories [12, IH, H, EE lill- From this viewpoint, 
the occurrence of a plateau structure at the early stage of nucleus-nucleus collision should not come as a surprise. 

The rapidity plateau distribution differs from the rapidity distribution of the bulk matter which is found to have a 
Gaussian shape [bI, [6^, [G^I . It is important to note that jets occur at an early stage of the nucleus-nucleus collisions, 
whereas the bulk medium properties are measured at the end-point of the nucleus- nucleus collision. A significant 
dynamical evolution must have occurred between the early beginning of the nucleus- nucleus collision and the end- 
point of the nucleus-nucleus collision. Therefore the early parton momentum distribution near the beginning stage of 
the nucleus- nucleus collision needs not be the same as the bulk matter distribution at the end-point of the nucleus- 
nucleus collision. One expects that starting with a non-isotropic plateau rapidity distribution that is much elongated 
in the longitudinal direction, a collision of two partons with large and opposing longitudinal momenta in adjacent 
spatial locations will redistribute the partons from the longitudinal direction towards the transverse directions, with 
a decrease in the longitudinal momenta of the colliding partons. Hence, the evolution will likely smooth out the 
anisotropic plateau rapidity structure to a significant degree as time proceeds. 

The subject of our focus being the near-side ridge and jet quenching in the early collision history, how the parton 
distribution function i^(r, p, t) evolves subsequently from the initial state to the end-point of nucleus-nucleus collision 
is beyond the scope of the present manuscript. The complete problem of parton evolution is a problem of great 
complexity [10, [U, [H, [13, HH, [H, [13 1 involving perturbative and non-perturbative elements. For example, in one of 
the descriptions using the Color-Glass-Condensate treatment of the initial conditions [ll|, it is not well-understood 
even within the Color-Glass-Condcnsatc community how the initial large rapidity correlations can evolve into to a 
thermal distribution in a short period of time of 1-2 fm/c or to a Gaussian rapidity distribution at the end-point of the 
nucleus-nucleus collision. Some recent advances suggest intrinsic color plasma instabilities that can lead to a breaking 
of the boost invariancc [3ft . ,80,. i81j, and other investigations suggest the bottom-up scenario involving gg — > ggg 
(s^ . [83| . The extracted early momentum distribution of a rapidity plateau obtained here serves to high-light the 
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important and unsolved issues of parton evolution that arc left outstanding by the present findings of this manuscript. 

The momentum loss of the jet parton and the geometry of the jet trajectory are other important aspects of the 
momentum kick model. The magnitude of the momentum kick imparted onto the medium parton has been found to 
be g = 1.0 GeV per jet-(medium parton) collision. This momentum gain by the kicked parton is clearly related to 
the momentum loss of the jet as a result of the jet-(medium parton) collisions. One obtains a good phenomenological 
description of the experimental data of the centrality dependence and coUisional energy dependence of Raa and 
the ridge yield. The extracted physical quantities furnish important, albeit approximate, empirical data for future 
investigations on the dynamics of parton production, parton evolution, and jet energy loss. The subject will come 
over and over again, each time with more and more accuracy and refinement, as we go through our course in physics. 

The successes of the simplifying model indicates that the momentum kick model contains promising degrees of 
freedom for the description of the gross features of the ridge phenomenon and jet quenching. There is however 
a limited range for the application of a completely analytical formulation, as many refinements and improvements 
necessitate additional degrees of freedom. Among other things, we envisage the need for a better description of 
the elementary jet-(medium parton) collision process, a better description of the dynamics of the medium, and the 
inclusion of effects of medium transverse collective and elliptic fiows that depend on the reaction plane orientations 
and medium spatial locations. There is the further complication for intermediate pt trigger particles that some of 
the trigger particles may arise not from the jets but from the medium [s^. A Monte Carlo implementation of the 
momentum kick model that will allow the inclusion of many refinements and improvements, and will therefore be of 
great interest. 

The author wishes to thank Drs. Fuqiang Wang, V. Cianciolo, Jiangyong Jia, Zhangbu Xu, and C. Nattrass for 
helpful discussions and communications. This research was supported in part by the Division of Nuclear Physics, U.S. 
Department of Energy, under Contract No. DE-AC05-00OR22725, managed by UT-Battelle, LLC. 
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